X-ray absorption signatures of the molecular environment in water and ice 



Wei Chen\ Xifan Wu^ and Roberto Car^'^'^l 

^Department of Physics, Princeton University, Princeton, NJ 08544, USA 
^Department of Chemistry, Princeton University, Princeton, NJ 08544, USA 
^Princeton Institute for the Science and Technology of Materials, 
Princeton University, Princeton, NJ 08544, USA 
(Dated: September 21, 2009) 

The x-ray absorption spectra of water and ice are calculated with a many-body approach for 
electron- hole excitations. The experimental features, including the small effects of temperature 
change in the liquid, are quantitatively reproduced from molecular configurations generated by ab- 
initio molecular dynamics. The spectral difference between the solid and the liquid is due to two 
major short range order effects. One, due to breaking of hydrogen bonds, enhances the pre-edge 
intensity in the liquid. The other, due to a non-bonded molecular fraction in the first coordination 
shell, affects the main spectral edge in the conversion of ice to water. This effect may not involve 
hydrogen bond breaking as shown by experiment in high-density amorphous ice. 
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The nature of the hydrogen bond (H-bond) network in 
water continues to be at the center of scientific debate [l| . 
A few years ago high-resolution x-ray absorption spectra 
(XAS) probed the local order of liquid and crystalline 
water phases, including the effect of temperature in the 
liquid [2[. These experiments stirred a storm of contro- 
versy because of their interpretation suggesting that a 
large fraction (>80%) of H-bonds are broken in the liq- 
uid [2I . This would imply an environment consisting pri- 
marily of chains and rings of H-bonds, in stark contrast 
with the conventional near-tetrahedral picture supported 
by diffraction [3] and other thermodynamic and spectro- 
scopic data [4, 5]. The broken H-bond fraction was esti- 
mated from the intensity of the pre-edge, which is promi- 
nent in the liquid and was believed to be absent in bulk 
ice. Subsequent x-ray Raman spectra, however, showed 
that the pre-edge feature is present, albeit with different 
intensity, not only in the liquid, but also in hexagonal 
(Ih), cubic (Ic), low-density amorphous (LDA) and high- 
density amorphous (HDA) ice [6]. This experiment re- 
ported another interesting observation, namely that wa- 
ter and HDA ice have spectra with the main-edge more 
prominent than the post-edge, while the opposite behav- 
ior occurs in Ih, Ic, and LDA ice. This is puzzling given 
that both LDA and HDA ice are disordered structures 
with an insignificant fraction of broken H-bonds. Spec- 
tral calculations, based on electronic density-functional 
theory (DFT) and near-tetrahedral liquid models, cor- 
rectly predicted the presence of three spectral features 
in ice and water and associated the enhancement of 
the pre-edge intensity to broken H-bonds 0, 0, 0, [13] • 
However, the agreement with experiment was only semi- 
quantitative and the calculations did not identify the 
cause of the significant changes in main- and post-edge 
spectra. These effects were generically attributed to dis- 
order 0, but this does not explain why LDA and HDA 
ice show opposite behavior 0. Finally, no attempt was 
made to discuss the effect of temperature in the liquid 



spectra, for which contrasting experimental data have 
been reported [i, [HI, H . 

In this paper we use liquid structures generated by 
ab-initio MD to compute x-ray spectra by adopting the 
final state rule [14] within a many-body formalism for 
electronic excitations [15]. Our calculations reproduce 
quantitatively the experimental spectra and their tem- 
perature dependence, supporting prior theoretical claims 
that the spectra are fully consistent with the conventional 
near-tetrahedral model of the liquid [8] . We associate the 
x-ray features to molecular excitons strongly influenced 
by the medium. Broken H-bonds do enhance the pre- 
edge intensity as pointed out in Ref. [2], but the differ- 
ence between water and ice is also due to other factors 
commonly associated to liquid disorder and does not re- 
quire an exceedingly large fraction of broken bonds 
A novel result of our study is the finding that the strong 
post-edge feature of ice originates from a peak in the den- 
sity of states that disappears in the liquid, as required to 
explain the corresponding spectral changes. This behav- 
ior is caused by the presence of a non-bonded fraction 
within the first molecular coordination shell, an effect not 
directly associated to broken H-bonds, even though the 
short-range order (SRO) of the liquid reflects the partial 
collapse of the H-bond network. We argue that the same 
effect, in absence of H-bond breaking, should explain the 
post-edge feature of HDA ice. 

Our approach can be viewed as an approximate version 
of the Bethe-Salpeter equation (BSE) for electron-hole 
excitations [l5| , which was recently applied with success 
to compute the optical absorption spectra of Ih ice [l^ 
and liquid water In the present case, we neglect 

the dynamics of the oxygen Is core hole, which is local- 
ized and long lived on the time scale of the absorption 
process. Then the BSE formalism reduces to a quasi- 
particle equation for the excited electron in presence of 
the frozen core hole and the sea of the other electrons [l3] • 
We solve this equation within the GW approximation, 



more specifically, the static Coulomb-hole and screened 
exchange (COHSEX) approximation [l5j], which is justi- 
fied for quasi-particle energies in the near edge range of 
the absorption spectrum. In this approach the local ex- 
change correlation potential of DFT is replaced by a non- 
local self-energy operator, which we evaluate by adopting 
a real space representation based on Maximally Localized 
Wannier functions H, [2^ . 



In Fig. [T] we compare calculated and experimental spec- 
tra Q of ice and water. The spectra are aligned at the 
onset and all the theoretical curves are multiplied by the 
scaling factor that adjusts the intensity of the main wa- 
ter peak (at T=363K) to the corresponding experimen- 
tal value (at STP). We use uniform Gaussian broadening 
with a full width at half maximum of 0.6 eV and 0.4 eV, 
for ice and water, respectively, to remove discreteness 
from the calculated spectra [2l| and to account for the 
ignored quasiparticle lifetime and intrinsic broadening ef- 
fects. The liquid spectra are extracted from a single snap- 
shot of 32 molecules at each temperature. We checked 
with calculations at the DFT level of theory that includ- 
ing more snapshots in the average has minor effect, in 
agreement with Ref. [sj] . At the two temperatures of the 
simulation the molecules are approximately tetrahedrally 
coordinated with ~ 7% and ~ 11% of broken H-bonds at 
T=330 K and T=363 K, respectively, according to the 
criteria of Ref. [22], or with - 14% and - 18% broken H- 
bonds, respectively, according to the criteria of Ref. [2]. 

The three main experimental features, pre-edge (533 
to 536 eV), near-edge (537 to 539 eV) and post-edge 
(from 539 eV and beyond), are well reproduced in terms 
of position, intensity, and spectral width. By contrast, 
previous DFT calculations, using excited- [sl or full-core 
hole potentials, underestimated the overall spectral 
width by ~2 eV or more, respectively, with correspond- 
ing shifts in the relative peak positions. Somewhat better 
spectral widths were obtained in half-core hole calcula- 
tions [10|, l23|, but this approach gave worse agreement 
in the pre- and main-edge 0, 0] • These inaccuracies are 
caused by the local DFT exchange-correlation potential 
that acts equally on all the excited states 0. By con- 
trast the non-local self-energy operator causes higher en- 
ergy states to experience reduced exchange effects, i.e. a 
reduced attractive potential, for an overall increase of the 
spectral width. The agreement between theory and ex- 
periment is very good in the liquid, in which we average, 
at each snapshot, over the 32 local configurations corre- 
sponding to the possible oxygen excitations in the super- 
cell. Use of a single high-symmetry local configuration 
in ice is responsible for the splitting of the post-edge fea- 
ture into two relatively sharp peaks, at lower and higher 
energy, respectively, compared to the single, broader, ex- 
perimental feature [lH. In the liquid (Fig. [T] panels (c) 
and (d)), a temperature increase enhances pre- and near- 
edge at the expense of post-edge intensity. Qualitatively, 
the same effect occurs in the transition from ice to wa- 
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FIG. 1: Calculated and experimental [2] XAS spectra. From 
bottom to top. (a) Ice: full line (theory, ice Ic), dashed line 
(experiment, ice Ih, at 95K). (b) Water: full line (theory 
at 363K), dashed line (experiment at 290K). (c) Theoreti- 
cal spectra for water at two different temperatures: blue line 
at 330K, red line at 363K. (d) Difference spectra: theory (full 
line), experiment (data points with error bars). The temper- 
ature difference in the experiment (AT = 65 K) was nearly 
double that in the calculation (AT = 33K). Consequently 
the theoretical data are magnified by a factor of 10 while the 
experimental data are magnified by a factor of 5 instead than 
by the factor of 10 used in Ref. [2]. A nearly linear dependence 
on temperature is to be expected for a temperature change 
small on the energy scale of the structural changes. 



ter [5]. Comparison of theoretical and experimental dif- 
ference spectra in panel (d) shows that the small effects 
of a temperature change are accurately reproduced, as il- 
lustrated by the isosbestic point that falls at 538.4 eV in 
the calculation, and at 538.8 eV in the experiment. Such 
a close correspondence in a derivative property would be 
unlikely barring a good correspondence of the trends in 
the molecular environment. A different behavior in which 
both main- and post-edge features decrease with temper- 
ature was reported in nanodroplet experiments |lll |. 

In our calculation, the intensity of the pre-edge fea- 
ture is weaker than experiment in both water and ice, 
but the intensity ratio is reasonably reproduced. This 
reffects approximations as our DFT calculations show 
that vibrational effects in ice and a larger simulation cell 
in water would both enhance the pre-edge. An overall 



underestimate of the pre-edge feature could also result 
from the assumption of a fully screened core- hole. More- 
over the network distortion and the fraction of broken 
H-bonds in a simulation using classical nuclei and the 
PBE exchange-correlation functional [26] might still be 
somewhat underestimated in s pite of a temperature ~60 
K higher than experiment 13, 27| . 

In Fig. [2] we report the COHSEX quasiparticle wave 
functions of representative pre-edge, near-edge and post- 
edge states of ice and water. We also report the three low- 
est unoccupied Hartree-Fock states of the water monomer 
with a fully relaxed oxygen Is core hole. The latter are 
antibonding states of 4ai, 2^2 and 862 symmetry, respec- 
tively, and correspond to the three lowest energy features 
of the XAS spectrum in the gas phase ^] . The similar- 
ity of condensed and gas phase orbitals suggests that the 
XAS features in ice and water should be attributed to 
molecular excitons strongly perturbed by the environ- 
ment. 



m-1 




W1 % 


m-2 




w-2 ^ 


m-3 







FIG. 2: Excited electron wave function in presence of a fully 
relaxed core hole on the central oxygen atom. The calcula- 
tion is at the Hartree-Fock level for the monomer (m), at the 
COHSEX level for ice Ic (i) and water (w). 1: pre-edge, 2: 
near-edge, 3: post-edge. Only molecules up to the second 
coordination shell are shown for ice and water. The color dif- 
ference in the density plots corresponds to a difference in the 
sign of the wave function. 

We assign the pre-edge to a bound exciton level whose 
electron orbital is similar to that of the valence exciton 
that marks the onset of optical transitions in ice [igI]. 
Although localized on the excited molecule, this orbital 
has significant weight also on the four nearest molecules. 



but the degree of localization may depend on our ap- 
proximations, such as e.g. fully screened core hole po- 
tential and partial neglect of local field effects in the 
COHSEX calculation. The intact H-bond network of ice 
enhances the symmetry and reduces the p-character of 
the excited quasi-electron on oxygen, lowering the dipo- 
lar transition amplitude and the pre-edge intensity, but 
the effect is not sufficient to make the feature disappear. 
This would require the distinction between donor and ac- 
ceptor bonds to vanish as in ice X, a stable form of ice 
under extreme pressure [29]. From detailed analysis we 
find that, while on average configurations with broken 
H-bonds contribute more to the pre-edge intensity than 
configurations with intact bonds, some distorted tetra- 
hedral configurations can be equally effective. Distortion 
and fraction of broken H-bonds increase with tempera- 
ture, enhancing the pre-edge intensity. 

The remaining two features in the near- and post-edge, 
correspond to more delocalized states and originate from 
several closely spaced levels in our condensed-phase su- 
per cell calculations, consistent with the assignment to 
exciton resonances. While in the monomer and the liq- 
uid higher energy states are less localized, the opposite 
happens in ice, as also reflected in the corresponding rel- 
ative intensities of near- and post-edge peaks. This phe- 
nomenon is not caused by the transition matrix element 
as both near- and post-edge resonances have strong oxy- 
gen p-character, but originates from the density of states 
(DOS). This is illustrated in Fig. [3] which depicts the DOS 
of water and ice (in absence of core hole within DFT). 
The conduction DOS of ice has a sharp peak at ~8 eV 
above the band edge corresponding to antibonding molec- 
ular orbitals of 62 symmetry. These orbitals have lobes 
protruding from the covalently bonded hydrogens and on 
the lone pair side of oxygen. In ice the antibonding state 
of one molecule does not overlap with the corresponding 
state of an adjacent molecule, originating a sharp peak 
in the DOS (Fig. [3]). Due to the distinction between an- 
tibonding and bonding directions, the overlap between 
antibonding states of H-bonded molecules is effectively 
forbidden. By contrast, when a non-bonded molecule is 
present in the flrst molecular coordination shell, overlap 
occurs, often involving the antibonding H lobes of the 
two adjacent non-bonded molecules (Fig. [3]). A substan- 
tial structural change is necessary for that to happen. In 
the liquid it follows from the partial collapse of the H- 
bond network. Orbital interaction broadens substantially 
the peak in the DOS shifting its center of gravity to lower 
energy, as shown in Fig. [3l As a consequence, post-edge 
intensity is reduced and oscillator strength is transferred 
to near-edge. This happens when ice melts and, to a 
lesser extent, when the temperature is raised in the liq- 
uid. Interestingly, the simulations of Ref. [8|] found that, 
when the liquid is quenched to form a glass, substan- 
tial post-edge intensity is recovered together with good 
tetrahedral coordination, albeit in a non-crystalline lat- 



tice. This structure should mimic LDA ice, which is ob- 
tained experimentahy by quenching the hquid. LDA ice 
has SRO close to cubic and hexagonal ice [30]. Not sur- 
prisingly, the corresponding x-ray spectra are very simi- 
lar [6]. On the other hand, HDA ice, which is obtained by 
pressure amorphization, sports a main edge similar to the 
liquid with a pre-edge closer to the crystal . The latter 
is consistent with the absence of broken bonds [s^ . Yet 
HDA has SRO different from crystalline ice, and actually 
closer to the liquid, due to the presence of an additional 
non-bonded molecule, on average, at non-tetrahedral lo- 
cations within the first coordination shell [30] . Arguably 
this is the reason for the similarity of the main edge of 
HDA ice and water, which both have molecular density 
higher than crystalline ice. In our liquid snapshots more 
than 10 percent of the molecules have overlapping an- 
tibonding orbitals. Interestingly, a sizeable fraction of 
these molecules (more than 60 percent) have five neigh- 
bors. Some of these local configurations correspond to a 
tetrahedrally bonded central molecule (2 donors, 2 accep- 
tors) with a non-bonded molecule additionally present in 
the coordination shell. Occasionally, the latter molecule 
may also be tetrahedrally bonded resulting in a local con- 
figuration similar to those of HDA ice [30] . 



tetrahedral picture. The calculated effects of a temper- 
ature change in the liquid are illuminating. Theory au- 
tomatically guarantees the relative normalization of the 
spectra at different temperatures, a condition that may 
not be straightforward to enforce in experimental spec- 



tra 



I2I . [l3l] . The excellent agreement of our calculation 
with the differential spectra of Ref. provides indepen- 
dent support of these data. It also suggests that, in spite 
of the shortcomings of current DFT approximations, the 
description of the H-bond network provided by ab-initio 
MD is fundamentally correct and can capture tempera- 
ture trends with surprising accuracy. Finally, we note 
that recent x-ray emission spectra (XES) of liquid water 
were interpreted as suggesting two coexisting local envi- 
ronments, one low coordinated and the other tetrahedral 
(sij . We are unable to directly address this issue because 
XES cross sections involve core-hole dynamics, which is 
beyond the present formalism, but we note that water 
structures from ab-initio MD do not seem to support a 
two-liquid hypothesis. 
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FIG. 3: (color online). Electronic DOS of ice (Ic) and wa- 
ter, calculated with a 4x4x4 A;-point mesh (ice) and at /c = 
(water) with uniform Gaussian broadening of 0.3 eV. The 
conduction band states have positive energy. The insets show 
a central tagged molecule and two adjacent molecules in ice 
(top) and water (bottom). Also plotted is the integrated lo- 
cal DOS in correspondence with the DOS peaks. The bottom 
inset shows overlap, indicated by an arrow, between the anti- 
bonding density of the non-bonded molecule on the left and 
the central molecule, while no overlap occurs between the lat- 
ter and the bonded molecule at its right. See also Ref. [lO] 

In conclusion, spectral calculations beyond ground- 
state DFT bear close correspondence to experiment. Our 
study confirms the sensitivity of the XAS spectra to SRO 
modifications due to structural transitions and temper- 
ature changes. H-bond breaking mainly affects the pre- 
edge while near- and post-edge are sensitive to an in- 
crease of coordination due to a non-bonded molecular 
fraction. Our work further confirms that the XAS spec- 
trum of water is consistent with the conventional near- 
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